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ABSTRACT 

Context. This work is part of a systematic X-ray survey of the Taurus star forming complex with XMM-Newton. 
«j Aims. We study the time series of all X-ray sources associated with Taurus members, to statistically characterize their X-ray variability, and 
compare the results to those for pre-main sequence stars in the Orion Nebula Cluster and to expectations arising from a model where all the 
X-ray emission is the result of a large number of stochastically occurring flares. 

Methods. The analysis of the lightcurves is based on a maximum likelihood algorithm that segments the time series in intervals of constant 
signal without the need of binning. Flares are defined with criteria that take into account the amplitude and the derivative of the segmented 
lightcurves. Variability statistics are evaluated for different classes of pre-main sequence stars (protostars, cTTS, wTTS, brown dwarfs), and for 
Q ' different spectral type ranges. Flare frequency and energy distribution are computed. 

^ Results. We find that roughly half of the detected X-ray sources show variability above our sensitivity limit, and in ~ 26 % of the cases this 
CZ2 variability is recognized as flares. Variability is more frequently detected at hard than at soft energies. The variability statistics of cTTS and 
wTTS are undistinguishable, suggesting a common (coronal) origin for their X-ray emission. The frequency of large flares (E > 10^^ erg) on 
^ , Taurus members is 1 event per star in 800 ksec. The typical duration of these flares - probably biased by the finite observing time - is about 
T ' 10 ksec. We have for the first time applied a rigorous maximum likelihood method in the analysis of the number distribution of flare energies on 
K> ' pre-main sequence stars. In its differential form this distribution follows a power-law with index a = 2.4 ± 0.5, in the range typically observed 
on late-type stars and the Sun. 

- - ' Conclusions. The signature of the X-ray variability in the pre-main sequence stars in Taurus and Orion provides twofold support for coronal 
heating by flares: (i) The correlation between the maximum variability amplitude and the minimum emission level indicates that both flare and 
quiescent emission are closely related to the coronal heating process, (ii) The power-law index a derived for the flare energy distribution is 
large enough to explain the heating of stellar coronae by nano-flares (a > 2), albeit associated with a rather large uncertainty that leaves some 
doubt on this conclusion. 
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1. Introduction 

X-ray variability is a characteristic property of magnetically ac- 
tive stars, including such diverse objects as the Sun, dMe flare 
stars, RSC Vn variables, and pre-main-sequence (pre-MS) T 
Tauri stars (Feigelson & Montmerle L 'l999: Favata & Micela L 
12003: Gtidel , 2004) . Variability on long time-scales may com- 
prise modulation by rotating active regions (with time-scales 
on the order of days corresponding to the rotation period of the 
star) or even activit y cycles (with expected time-scales on the 
order of years); see lMarino et al.l (l2003h : Ipavata et al. I (l2004l) 
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for examples. Most of the observed variations have short time- 
scales (on the order of hours), are accompanied by an increase 
in electron temperature, and are therefore attributed to flares 
resulting from magnetic reconnection events. 

The intense X-ray radiation from flares is ascribed to the 
filling of magnetic loops with heated gas that has been driven 
from the chromospheric layers along the magnetic field lines 
into the corona (the 'chr omospheric evaporation scenario', e.g., 
lAntonucci et al.Lll984l) . Various authors have therefore spec- 
ulated that the overall X-ray emission seen in magnetically ac- 
tive stars is in fact the result of a large number of stochastically 
occurring flares. In analogy t o the 'mic roflare heating' hypoth- 
esis in solar physics (see e.g. lHuds'on ), these events are 
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believed to be distributed in energy such that the large majority 
releases only small amounts of energy, making their individual 
identification in light curves impossible while we only measure 
the in t egrate d emission (e.g., Audard et al. , 2000; Giidel et al. , 
l2002l[^)()l . 

Extensive studies yielding statistical properties of X-ray 
variability for different classes of active stars have long been 
impeded, mostly due to observational restrictions deriving from 
insufficient monitoring time. Further limitations, besides sen- 
sitivity, are imposed by the poor temporal coverage of X-ray 
instruments in low-Earth orbits that result in frequent occul- 
tations of the target by the satellite revolution. XMM-Newton 
and Chandra have opened the field for systematic investiga- 
tions on time-scales typical for magnetic activity. This is due 
to the eccentric orbits of both satellites, enabling uninterrupted 
observations of up to ~ 2 d. Even then, variability studies are 
difficult, because in practice it is nearly impossible to maximize 
both sample size and exposure time per target. An outstanding 
study in this respect is the nearly 13 d long observation of the 
Orion Nebula Cluster (ONC) by Chandra, termed the Chandra 
Orion Ultradeep Project (COUP). The COUP yielded not only 
the longest uninterrupted X-ray time series for pre-MS stars so 
far, but owing to the dense population of the ONC it allowed 
also to examine in a homogeneous way the X-ray variability 
of the largest sample of pre-MS stars: More than 1600 X-ray 
sources were detected in a single Chandra field, mostly mem- 
bers of the ONC. Amongst others, the COUP resulted in a sys- 
tematic study of flares on 'young Suns ' (Wolk et al. , 2005), 
rotational mo dulation i Flaccom io et aTl EoOSi) . the geometry 
of flare loops (^Fa vata et al. i.i2005l) . and X-rav variability of hot 
stars ( Stelzer et al. L I2005I) . 

In contrast to the ONC, the Taurus st ar forming c omplex 
spans a large angle on the sky (see Fig. 1 in ICiidel et al.[l2006) . 
This implies that observing a significant fraction of its stellar 
members with a telescope operating in pointing mode requires 
many different exposures. Realistically, the whole population 
can b e sampled only by an all-sky monitoring. Neuhauser et al. 

have used the data obtained during the ROSAT All-Sky 
Survey (RASS) to study the X-ray emission of the pre-MS 
stars in Taurus. However, the RASS consisted in snapshots of 
a few seconds each separated by ~ 90min, corresponding to 
the Earth revolution of the ROSAT satellite, and is inappropri- 
ate for timing studies. A systematic study of the X-ray vari- 
ability of the pre-MS stars in Taurus-Auriga based on more 
than 100 pointed RO SAT PSPC observations was presented by 
IStelzer et al. 1 (I2OOOI) . That investigation yielded the first esti- 
mate of flare rates for these stars, which were found in outburst 
during about 1 % of the observing time. Various uncertainties 
compromise this measurement, e.g. the ROSAT data set was (i) 
inhomogeneous in sensitivity due to the vastly different expo- 
sure times of the individual observations, (ii) affected by Earth 
blocks due to the satellite's revolution, and (iii) limited to soft 
energies. 

Here we present a variability study for the pre-MS stars 
in the Taurus star forming region based on extensive XMM- 
Newton observations. We discuss data from the XMM-Newton 
Extended Survey of the Taurus Molecular Cloud (XEST). This 
project was devised to comprise observations with roughly uni- 



form and uninterrupted exposures of ~ 30 - 40ksec. Thus, it 
represents a major improvement with respect to previous X-ray 
data of the Taurus Molecular Cloud (TMC). A further advan- 
tage is the extension of the hard energy band beyond 2keV. 
A detailed descrip tion of the XEST observations is found in 
iGiidel et'al] (I2OO6V Here we present a systematic time series 
analysis of the pre-MS population detected in the XEST. 

In Sect. 121 the methods used in the analysis of the time se- 
ries are described, and the results are summarized. Sect. ideals 
with a quantitative investigation of the lightcurves. The rela- 
tively short duration of the XEST observations imply that the 
most likely cause of the variability is flaring. We introduce our 
way of defining flares and present the results from the flare 
detection process in that section. Variability statistics for dif- 
ferent classes of pre-MS stars in the TMC are given in Sect. ^ 
including a discussion of observational biases. The frequency 
and the energy distributions of the detected flares are examined 
in Sect.lsjand Sect.|6l Sect.0presents a detailed comparison of 
the variability characteristics of cTTS and wTTS. We summa- 
rize our results in Sect.|Sl 

2. Time series analysis 

The analysis presented in this paper includes the TMC mem- 
bers of all XMM-Newton fields from the XEST, i.e. the obser- 
vations from our survey with XMM-Newton and the observa- 
tions added from the XMM-Newton archive (see observing log 
in Table 1 of Gudel et al., 2006), except field No. 1 (TTau) and 
No. 25 (AA Tau), that are dedicated to separate projects (Giidel 
et al., in prep.; Grosso et al., in prep.). Some fields are spa- 
tially overlapping, such that several stars are detected in more 
than one exposure. Especially, in XEST fields No. 23 and No. 24 
the satellite was pointing at the same sky position and the two 
exposures are adjacent in time (separated by only ~ Iksec). 
Therefore, their data was merged for the time series analysis. 
The total number of X-ray sources in all studied XEST fields 
that are identified with known TMC members after merging of 
XEST-23 and XEST-24 is 126. 

The time series analysis is based on the EPIC/pn photon 
event lists for the individual X-ray sources. There are three ex- 
ceptions where we resort to the MOS data: (i) for XEST field 
No. 26 because there is no EPIC/pn data, (ii) for sources that are 
outside the field-of-view of EPIC/pn, or (iii) sources located on 
a chip gap resulting in a major loss of photons in EPIC/pn. 

The photon extraction regions for source and background 
photons have been defined taking care to avoid contributions 
from adjace nt sources. The procedure is described in detail by 
iGiidel et'al] (I2006.) . For each X-ray source detected in both 
XEST-23 and XEST-24 the event lists of the source and of 
the background of both observations were merged. For a given 
source the original photon extraction areas in XEST-23 and in 
XEST-24 may be slightly different. To avoid introducing artifi- 
cial variability we selected for each star the smaller of the two 
extraction regions for both exposures. 

The statistical methods we used are independent of data 
binning, and are described below. The analysis was carried out 
in different energy bands: 0.3 - 7.8 keV (broad), 0.3 - l.OkeV 
(soft), and 1.0 - 7.8 keV (hard). 
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2.1. Maximum likelihood blocks 

We searched for variability in each photon time series using a 
ma ximum hkehhood a lgorithm. The technique was described 
by IWolk et al~l ( l2005l) . who have used the same method for 
an X-ray variability study of a subsample of the ONC mem- 
bers detected in the COUP. It is derived from the Bayesian 
blocks described by Scargle (1998). In brief, the time series 
is split into periods of constant signal (so-called 'segments' or 
'blocks') under the assumption of Poisson noise, by search- 
ing iteratively for change points in the intensity level. In con- 
trast to conventional lightcurve analysis this method works di- 
rectly on the sequence of photon arrival times without the need 
for binning, such that ambiguities introduced by the choice 
of bin size and bin start are avoided. The maximum likeli- 
hood block (MLB) algorithm has two free parameters, the min- 
imum number of counts per segment (A^min) and the confi- 
dence level (CL) for the established change points. The sig- 
nificance thresholds for the change points were determined 
for any given A^^in through extensive simulations of constant 
lightcurves (Flaccomio et al., in prep.). The results presented 
throughout this paper are for CL - 99.9% and A^min - 20. 
Now we proceed to a justification of this choice and a detailed 
description of the data analysis steps. 

While the background of Chandra is very low and can be 
approximated by a constant, in XMM-Newton observations the 
background is non-negligible and variable. This complicates 
the analysis. To take account of the background and its vari- 
ability we have extended the method with respect to the use for 
the COUP. In short, the principle of this procedure is to create 
a background subtracted 'source-only' or 'net source' events 
list, by removing individual photons from the events list of the 
source position. Subsequently, the variability tests can be ap- 
plied to this net source events Ust. 

In practice, for the treatment of the background we pro- 
ceed for each source events list and the associated background 
events list in the following way: First, the MLB algorithm is 
applied to the time series for the background, i.e. the back- 
ground events list ('B') is divided into ML blocks by identi- 
fying change points in the background level. The matter is then 
how to subtract the background from the events detected at 
the source position (which at this point comprise both source 
and background photons; - 'S-hB"). We first scale the seg- 
mented background to the extraction area of the source. This 
gives the expected number of background events at the source 
position for each segment in which the background is con- 
stant. The respective number of photons is then subtracted from 
the 'S-hB" events list uniformly within each time-interval of 
constant background. The result is a background subtracted 
source-only events list ('S'), ready for further analysis. 

The MLB algorithm is now applied to 'S' to yield the seg- 
mented source lightcurve. When computing the errors of the 
count rate in each segment, the (subtracted) number of back- 
ground photons within the respective time segment is consid- 
ered. Note that, obviously, the segments of 'S' are different and 
independent from the segments of 'B'. 

In general, a source can be considered variable at a cer- 
tain confidence level if the background-subtracted blocked 



lightcurve consists of more than one segment. In the course 
of the analysis it was noticed that for some sources the de- 
tected variability coincides with times of very high and vari- 
able background. These variations could be a result of inac- 
curate background subtraction. To check if the variations in 
these critical time intervals are spurious we simulated constant 
lightcurves by generating a random distribution of photon ar- 
rival times within a time interval corresponding to the length 
of the observation. The number of photons randomly gener- 
ated was chosen such that the count rate of the simulated data 
is equal to the average net count rate of the observed source. 
Then, we superposed the observed area-scaled background on 
the simulated source events list, and carried out the background 
subtraction and MLB analysis. The procedure for the treatment 
of the background is analogous to the one described above for 
the observed data, i.e. the background was subtracted block- 
wise from the simulated 'Ssim + B' events list and a background 
subtracted simulated 'Ssim' events list was computed. 

For each source 5000 simulations were performed, with 
CL = 99.9 % and Nrmn = 20. Then the blocked lightcurves 
were examined. The simulated data sets are expected to yield 
no change point in the MLB analysis. However, in a substantial 
number of sources where the algorithm detected variability in 
the observed data, variability was also detected in a high frac- 
tion (» 10%) of the simulated lightcurves. These cases are 
found among faint sources that have high and strongly variable 
background, and the statistical subtraction of background pho- 
tons is problematic. 

To avoid such detections of spurious variability, for the ul- 
timate analysis of the data we excluded all time intervals cor- 
responding to blocks from 'B' in which the number of events 
in the background after scaling to the source area is three or 
more times higher than the number of events in the source. 
In general, the reduction of the exposure time by this 'high- 
background filter' is low; less than 10 % for 105 of 126 sources. 
The worst cases occur in XEST-03 and XEST-24, where the 
background is so high, that - depending on the source bright- 
ness - up to 40 % of the observing time was removed. 

Residual variability imposed by the background can be as- 
sessed with the help of the simulated data. Clearly, the effect of 
background induced variations is to decrease the significance 
of any detected variability with respect to the chosen CL. In 
practice, we use the fraction of simulated event lists that ap- 
pear non-variable with the MLB method to define a 'corrected' 
confidence level CLcon for the variability detected in the ob- 
served data. We delay the discussion of the results to Sect. ^ 
and proceed first with the description of the analysis steps. 

The analysis was performed in all three energy bands (soft, 
hard, and broad) with values for A^min of 1 and 20, henceforth 
referred to as MLB 1 and MLB 20 respectively. Higher values 
for Nmin can be used e.g. to define time segments for time re- 
solved spectral analysis (see Franciosini et al. , 2006). A low 
number of the allowed minimum counts per block favors the 
detection of variability in faint sources, but also the detection 
of spurious variations, that arise e.g. from residual background 
contamination. A lightcurve with an example for this latter ef- 
fect is shown in Fig.[n see figure caption for further explana- 
tions. 
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Fig. 1. Broad band EPIC/MOS 1 lightcurve of XEST-19-009. The background subtracted source signal is shown in two represen- 
tations: binned into 1000 s intervals (red) and segmented with the MLB algorithm (blue). The left panel shows the segmentation 
obtained with A^min - 20 (only one block, i.e. constant lightcurve), the right panel shows the segmentation resulting from A^min - 1 
(two blocks, i.e. variable lightcurve). The background segments scaled to the source area are shown in yellow. Time intervals 
rejected by our high-background filter are represented by grey-shaded areas. The division into two blocks with MLB 1 probably 
results from imperfect background subtraction. 



A careful comparison showed that the variability statistics 
obtained with MLB 1 and with MLB 20 are widely consistent. 
For seven sources the segmentation of the broad band is dif- 
ferent for MLB 1 and MLB 20. Among these, three sources are 
constant for A^min = 20, but one change point is found with 
MLB 1 . Since these three sources are very faint, the detection 
of variability is hampered with MLB 20 (because the total num- 
ber of net source counts is only slightly higher than the change 
point threshold A^min - 20). However, in practice, their variabil- 
ity detected with MLB 1 is likely spurious because related to 
times of high background not rejected by our high-background 
filter; see example in Fig.^ For the other four sources at least 
one change point is detected with both MLB 1 and MLB 20, 
and only the number or position of the change points is differ- 
ent between the two tests. To conclude, we find no case where 
the detection of variability is impeded with MLB 20 due to the 
higher threshold for the minimum number of counts per seg- 
ment with respect to MLB 1. Based on these results, we de- 
cided to present in this paper the results of the analysis with 



2.2. Kolmogorov-Smirnov test 

As an additional investigation of the variability we performed a 
Kolmogorov-Smirnov (KS) test on the background subtracted 
events lists. The output of the KS test is a probability Pks that 
the variability detected in the observed data is physical. Since 
the 'S' events file has been manipulated by the subtraction of 
the background, it is not immediately obvious that the KS-test 
is applicable. To examine the validity of the KS analysis we 
made use of the simulations of constant source event lists de- 
scribed above. Analogously to the MLB analysis, the simulated 
lightcurves should be non-variable against the KS-test, except 



for spurious variability detections due to random fluctuations in 
the photon arrival times and additional noise introduced by the 
background subtraction. To quantify the effect of such spurious 
variations, we computed for each variable source the fraction fp 
of simulated lightcurves that show a probability Pks of being 
variable higher than the value derived for the observed photon 
time series. The 'effective' probability for variability in the data 

is ^'KS.eff = 1 - /p- 

2.3. Results from variability tests 

The variability analysis described in the previous section was 
applied to 126 photon time series from the 122 different TMC 
members detected in the XEST. The results are given in Tabled 
for all examined sources. For a few sources, in the soft and/or 
hard band the variability tests could not be carrie d out due to 
poor statistics. The entries in Tableware given as in lGiidel et alJ 
( 2006) in order of increasing right ascension. We list the XEST 
No. (col. 1) , identification (col. 2), classification defined by 
iGiidel et ail (E0O6.) that combines the Young Stellar Object 
class derived from the infrared spectral energy distribution and 
the T Tauri type based on Ha equivalent widths (col. 3), and 
the spectral type (col. 4). More details concerning the defini- 
tions of the stellar parameters are found in Gudel et al] ll2006l) . 
Col. 5 gives the broad band exposure time. Generally, the in- 
strument used is EPIC/pn; the exceptions are indicated by la- 
bels 'Ml' for MOS 1 and 'M2' for MOS 2 at the end of col. 5. 
The remaining columns represent the results of the variability 
tests. For all three energy bands the number of source photons 
('Cts'), the numberof segments resulting from the MLB analy- 
sis ('A^b')i the confidence level of the MLB analysis ('CLcorr'). 
and the probability for variability according to the KS test 
(^KS.eff) are given; see table caption for more details. 



Table 1. X-ray variability of TMC members during the XEST. Cols. 1-4 give the XEST number, stellar identification, young star class, and spectral type (see iGiidel et'all 
2006, for details on these classifications). The summed broad band GTI after removal of time intervals with high background is given in col. 5. The analysis is based on 
EPIC/pn data; exceptions are marked with labels Ml for MOS 1 and M2 for MOS 2. The remaining columns represent the number of source photons ('Cts'), the number of 
segments i'Nb), the confidence level of the MLB test ('CLcoit'). and the probability for variable source according to the KS-test ('/"KS.eft')- Results are given for the broad 
(0.3 - 7.8keV), soft (0.3 - 1.0 keV) and hai'd (1.0 - 7.8keV) band. 
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2 


0.67 


1.00 


17 


1 




0.75 


178 


2 


0.74 


1.00 


23-037/24-032 


V410 X2 


9 


MO 


103309 


3625 


4 


0.73 


1.00 


475 


2 


0.78 


0.77 


3133 


3 


0.41 


1.00 


23-045/24-038 


V410 X4 


9 


M4 


94996 


810 


3 


0.60 


1.00 


89 


1 




0.71 


749 


3 


0.73 


1.00 


23-047/24-040 


V892 Tau 


5 


B9 


109954 


28548 


11 


0.97 


1.00 


2354 


2 


0.96 


1.00 


26196 


13 


0.97 


1.00 


23-048 


LR 1 


9 


K4.5 


43482 


68 


2 


0.92 


1.00 


- 


- 


- 


- 


83 


1 




0.99 


23-050/24-042 


V410 X7 


3 


M0.75 


110155 


3259 


6 


0.91 


1.00 


413 


3 


0.91 


1.00 


2846 


6 


0.93 


1.00 


23-056/24-047 


Hubble 4 


3 


K7 


109845 


30574 


4 


0.96 


1.00 


14776 


1 




0.66 


15791 


7 


0.95 


1.00 


23-061/24-054 


V410 X6 


3 


M5.5 


97683 


436 


2 


0.84 


1.00 


222 


2 


0.88 


1.00 


207 


2 


0.82 


1.00 


23-063/24-055 


V410 X5 


3 


M5.5 


104095 


2058 


3 


0.88 


1.00 


631 


2 


0.60 


1.00 


1439 


3 


0.76 


1.00 


23-067/24-058 


FQ Tau AB 


2 


M3/M3.5 


95824 


75 


1 




0.79 


39 


1 




0.90 


37 


1 




0.35 


28-100 


BP Tau 


2 


K7 


129661 


28835 


12 


0.91 


1.00 


15671 


8 


0.97 


1.00 


13126 


12 


0.93 


1.00 


23-074/24-061 


V819TauAB 


3 


K7 


106945 


9857 


4 


0.85 


1.00 


6132 


3 


0.95 


1.00 


3699 


6 


0.70 


1.00 


11-023 


2MJ042 13459 


3 


M5.5 


31441 


77 


1 




0.94 


53 


1 




0.85 


26 


1 




0.04 


11-037 


CFHT-Tau 10 


3 


M5.75 


28006 


10 


1 




0.59 


10 


1 




0.46 










21-038 


RYTau 


2 


Kl 


45030 


6864 


3 


0.91 


1.00 


991 


2 


0.87 


1.00 


5877 


2 


0.92 


1.00 


21-039 


HD 283572 


3 


G5 


44346 


94102 


4 


0.98 


1.00 


55032 


2 


0.98 


1.00 


39073 


4 


0.97 


1.00 



Table 1. continued 



Broad band Soft band Hard band 



XEST No. 


Identification 


Type 


SpT 


Expo [s] 




Cts 


Nb 






Cts 


Afb 


CLcOTT 




Cts 


A'b 




^KS.eff 


11-054 


Haro 6-5 R 


1 


K5 


35463 




64 


1 




0.69 




1 




0.96 


59 


1 




0.77 


11-057 


FS Tau AC 


2 


M0/M3.5 


37353 




2845 


5 


0.97 


1.00 


120 


1 




0.28 


2729 


5 


0.97 


1.00 


21-044 


LkCa 21 


3 


M3 


45474 




3539 


2 


0.88 


1.00 


2467 


1 




0.62 


892 


2 


0.80 


1.00 


1 1-079 


CFHT-Tau 21 


2 


Ml. 25 


34673 




164 


2 


0.92 


1.00 


9 


1 




0.88 


141 


2 


0.95 


0.99 


02-013 


FV Tau AB 


2 


K5/K6 


38702 




310 


1 




0.99 


7 


1 




0.90 


300 


1 




0.99 


02-016 


KPNO-Tau 13 


3 


M5 


38190 




279 


1 




0.78 


122 


1 




0.68 


153 


1 




0.99 


02-022 


DG Tau A 


2 


K6 


39027 




681 


3 


0.94 


1.00 


267 


1 




0.78 


384 


3 


0.95 


1.00 


15-020 


JH 507 


3 


M4 


27652 




896 


2 


0.91 


0.99 


632 


2 


0.96 


1.00 


260 


1 




0.98 


13-004 


GV Tau AB 


1 


K3-7 


26617 




296 


3 


0.88 


1.00 


50 


1 




0.42 


248 


3 


0.90 


1.00 


15-040 


DH Tau AB 


2 


Ml 


30891 




16216 


5 


0.98 


1.00 


7840 


4 


0.98 


1.00 


8375 


5 


0.97 


1.00 


15-042 


DI Tau AB 


3 


MO 


30891 




2049 


1 




0.98 


1340 


1 




0.60 


693 


2 


0.82 


1.00 


15-044 


KPNO-Tau 5 


4 


M7.5 


18419 




40 


1 




0.90 


45 


1 




0.70 


22 


1 




0.93 


14-006 


TO Tau A 

1 CI d i V 


2 


M0.5 


30377 




1079 


3 


0.94 


1.00 


179 


2 


0.95 


1.00 


903 


3 


0.95 


1.00 


13-035 


FX Tau AB 


2 


Ml 


30547 


M2 


162 


1 




0.94 


81 


1 




0.90 


81 


1 




0.89 


14-057 


DK Tau AB 


2 


K7 


31879 




2358 


1 




0.28 


1021 


1 




0.67 


1329 


1 




0.05 


22-013 


MHO 9 


3 


M4.25 


52864 




355 


1 




0.98 


287 


1 




0.95 


77 


1 




0.80 


22-021 


MHO 4 


4 


M7.1 


53756 




313 


1 




0.53 


209 


1 




0.86 


105 


1 




0.54 


22-040 


L1551 IRS5 

Xj X »^ ^ X XX vw^ ^ 


1 




51467 




84 


1 




0.77 


31 


1 




0.05 


57 


1 




0.85 


22-042 


LkHa 358 


2 


M5.5 


52730 




70 


2 


0.63 


1.00 










81 


2 


0.56 


0.84 


22-043 


HL Tau 


1 


K5 


54434 




1416 


2 


0.97 


0.99 


5 


1 




0.88 


1409 


2 


0.98 


0.99 


22-047 


XZ Tau AB 


2 


M2/M3.5 


54434 




8181 


5 


0.97 


1.00 


2894 


4 


0.98 


1.00 


5286 


6 


0.97 


1.00 


22-056 


L1551 NE 


1 




52096 




16 


1 




0.72 










20 


1 




0.72 


03-005 


HK Tau AB 


2 


M0.5/M2 


19581 




45 


2 


0.92 


0.98 


14 


1 




0.99 


27 


1 




0.71 


22-070 


V710 Tau BA 


2 


M0.5/M2 


54346 




3191 


1 




0.82 


1611 


1 




0.79 


1569 


1 




0.19 


19-009 


JH 665 


3 


M5.5 


34356 


Ml 


23 


1 




0.99 


7 


1 




0.60 


19 


1 




1.00 


22-089 


L1551 51 


3 


K7 


54443 




4739 


4 


0.97 


1.00 


3127 


3 


0.98 


1.00 


1603 


3 


0.95 


1.00 


22-097 


V827 Tau 


3 


K7 


56844 


MI 


5255 


2 


0.98 


1.00 


2593 


2 


0.98 


1.00 


2661 


2 


0.98 


1.00 


03-016 


Haro 6-13 


2 


MO 


23372 




203 


2 


0.86 


0.99 


27 


1 




0.28 


188 


2 


0.88 


1.00 


22-100 


V826 Tau 


3 


K7 


54443 




18357 


3 


0.98 


1.00 


11587 


2 


0.98 


1.00 


6767 


2 


0.97 


0.99 


22-101 


MHO 5 


2 


M6 


52251 




533 


1 




0.80 


418 


1 




0.46 


124 


1 




0.21 






J 








104 


I 






XI 

1 


\ 




n 61 


in 


I 






03-019 


V928 Tau AB 


3 


M0.5 


32024 




1545 


1 




0.97 


826 


1 




0.98 


721 


1 




0.80 


03-022 


FY Tau 


2 


K5 


33093 




1271 


4 


0.96 


0.99 


362 


1 




0.89 


910 


3 


0.96 


1.00 


03-023 


FZTau 


2 


MO 


34957 


Ml 


117 


1 




0.99 


54 


1 




0.95 


63 


1 




0.97 


17-002 


IRAS 04295+2251 


1 




27556 




149 


2 


0.95 


1.00 


7 


1 




0.97 


140 


2 


0.95 


1.00 


19-049 


UZTauE-HW(AB) 


2 


Ml/2/3 


32225 




1292 


2 


0.92 


0.91 


599 


1 




0.44 


704 


2 


0.91 


0.77 


17-009 


JH112 


2 


K6 


28498 




453 


1 




0.99 


143 


1 




0.62 


307 


1 




0.99 


03-031 


CFHT-Tau 5 


4 


M7.5 


32755 


Ml 


37 


1 




0.83 


6 


1 




0.54 


28 


1 




0.96 


04-003 


CFHT-Tau 5 


4 


M7.5 


29570 




124 


1 




0.72 


25 


1 




0.30 


99 


1 




0.92 


03-035 


MHO 8 


3 


M6 


21097 




78 


1 




0.99 


73 


1 




0.81 


41 


2 


0.93 


1.00 



Table 1. continued 



Broad band Soft band Hard band 



Abil JNo. 


Identification 


Type 


apl 


Expo [s] 


Cts 


A'b 


^ ^corr 


^KS.eff 


Cts iVb 




^KS.eff 


Cts 




^ ^corr 


^KS.eff 


04-009 


MHOS 


3 


M6 


29561 


73 


1 




0.63 


54 1 




0.78 


20 


1 




0.18 


04-010 


GH Tau AB 


2 


M1.5/M2 


33077 


Ml 42 


1 




0.96 


10 1 




0.78 


32 


1 




0.93 


04-012 


V807 Tau (SNab) 


2 


K7/M3 


30317 


1733 


3 


0.98 


1.00 


1207 1 




0.94 


525 


3 


0.98 


1.00 


18-004 


KPNO-Tau 14 


3 


M6 


27895 


481 


3 


0.94 


1.00 


62 1 




0.85 


418 


3 


0.96 


1.00 


04-016 


V830 Tau 


3 


K7 


30411 


9455 


5 


0.89 


1.00 


5618 4 


0.98 


1.00 


3839 


5 


0.83 


1.00 


17-027 


IRAS 04303-1-2240 


2 


M0.5 


29039 


2257 


3 


0.96 


0.95 


114 2 


0.92 


1.00 


2148 


1 




0.88 


04-034 


GITau 


2 


K7 


30814 


748 


1 




0.86 


269 1 




0.91 


480 


1 




0.62 


04-035 


OK Tau AB 


2 


K7 


30814 


1334 


2 


0.96 


1.00 


481 1 




0.79 


853 


2 


0.96 


1.00 


18-019 


IS Tau AB 


2 


K7/M4.5 


27970 


1000 


2 


0.97 


1.00 


318 1 




0.92 


683 


2 


0.97 


1.00 


17-058 


CI Tau 


2 


K7 


27575 


243 


2 


0.94 


1.00 


40 2 


0.98 


0.98 


204 


2 


0.95 


1.00 


18-030 


IT Tau AB 


2 


K2 


28666 


9304 


2 


0.98 


1.00 


1102 1 




0.56 


8194 


2 


0.98 


1.00 


17-066 


JH 108 


3 


Ml 


28847 


2539 


4 


0.97 


1.00 


993 3 


0.97 


1.00 


1546 


5 


0.95 


1.00 


17-068 


CFHT-BD Tau 1 


4 


M7.1 


26906 


86 


1 




0.99 


13 1 




0.81 


70 


1 




0.99 


09-010 


HO Tau AB 


2 


M0.5 


21904 


61 


1 




0.31 


52 1 




0.37 


19 


1 




0.82 


08-019 


FF Tau AB 


3 


K7 


34927 


1345 


3 


0.88 


1.00 


657 2 


0.94 


1.00 


680 


2 


0.95 


1.00 


12-040 


DNTau 


2 


MO 


29424 


5110 


3 


0.99 


1.00 


3072 1 




0.97 


2040 


2 


0.98 


0.99 


12-059 


CoKu Tau 3 AB 


3 


Ml 


29426 


10548 


1 




0.98 


3547 1 




0.56 


7000 


2 


0.98 


0.98 


09-022 


KPNO-Tau 8 


3 


M5.75 


29425 


1609 


3 


0.96 


0.99 


900 1 




0.96 


721 


1 




0.97 


08-037 


HQ Tau AB 


3 




41566 


Ml 3194 


2 


0.93 


1.00 


937 1 




0.99 


2255 


2 


0.91 


1.00 


09-026 


HQ Tau AB 


3 




29826 


4840 


4 


0.98 


1.00 


1519 2 


0.98 


1.00 


3320 


2 


0.98 


1.00 


08-043 


KPNO-Tau 15 


3 


M2.75 


35322 


6014 


1 




0.97 


1636 1 




0.74 


4384 


2 


0.89 


1.00 


09-031 


KPNO-Tau 15 


3 


M2.75 


26257 


276 


1 




0.42 


140 1 




0.63 


136 


1 




0.84 


08-048 


HP Tau AB 


2 


K3 


37574 


2488 


2 


0.82 


1.00 


688 1 




0.97 


1801 


2 


0.83 


1.00 


08-05 la 


HPTau/G3 AB 


3 


K7 


35141 


992 


1 




0.87 


451 1 




0.38 


556 


1 




0.27 


08-051 


HP Tau/G2 


3 


GO 


35574 


15377 


2 


0.98 


1.00 


5871 2 


0.98 


1.00 


9510 


2 


0.99 


1.00 


08-058 


Haro 6-28 AB 


2 


M2/M3.5 


35457 


306 


2 


0.88 


0.19 


144 1 




0.17 


160 


1 




0.59 


08-080 


CFHT-BD Tau 3 


4 


M7.75 


35678 


29 


1 




0.99 


5 1 




0.81 


19 


1 




0.99 


05-005 


CFHT-Tau 6 


4 


M7.25 


25463 


68 


2 


0.87 


0.97 


34 1 




0.43 


33 


1 




0.99 


05-013 


GN Tau AB 


2 


M2.5 


24825 


105 


1 




0.93 


13 1 




0.20 


93 


1 




0.98 


05-017 


IRAS 04365+2535 


1 




24037 


48 


1 




0.34 


— — 


— 


— 


46 


1 




0.30 


05-024 


IRAS 04369-H2539 


2 


K4 


28083 


296 


3 


0.85 


1.00 


11 1 




0.94 


286 


2 


0.83 


1.00 


07-011 


JH223 


3 


M2 


29922 


227 


1 




0.69 


146 1 




0.53 


79 


1 




0.52 


07-022 


Haro 6-32 


3 


M5 


29850 


157 


1 




0.75 


95 1 




0.16 


62 


1 




0.95 


07-041 


IRAS 04385-t-2550AB 


2 


MO 


28467 


185 


3 


0.95 


0.99 


12 1 




0.93 


175 


3 


0.96 


0.99 


10-017 


CoKuLk332/G2 AB 


3 


MO.5/2.5 


29960 


3753 


1 




0.44 


1381 1 




0.33 


2374 


1 




0.56 


10-018 


CoKuLk332/Gl AB 


3 


K7/M1 


30036 


420 


1 




0.52 


152 1 




0.67 


268 


1 




0.30 


10-020 


V955 Tau AB 


2 


K5/M1 


31665 


Ml 247 


1 




0.99 


28 1 




0.74 


219 


1 




0.98 


10-034 


CIDA7 


2 


M4.75 


25405 


54 


1 




0.83 


40 1 




0.24 


19 


1 




0.67 


10-045 


DP Tau 


2 


M0.5 


29926 


56 


1 




0.98 


18 1 




0.99 


34 


1 




0.99 


10-060 


GO Tau 


2 


MO 


28564 


256 


1 




0.78 


95 1 




0.75 


166 


1 




0.60 



Table 1. continued 



Broad band Soft band Hard band 



XEST No. 


Identification 


Type 


SpT 


Expo [s] 




Cts 


A^b 


^^corr 




Cts 


A^b 


^-^corr 




Cts Afb 


^^corr 


^KS.eff 


26-012 


2M J04552333+30 


4 


M6.25 


126315 


Ml 


25 


1 




0.65 


13 


1 




0.94 


12 1 




0.57 


26-034 


2M J04554046-I-30 


3 


M5.25 


126497 


Ml 


19 


1 




0.30 


10 


1 




0.02 


10 1 




0.91 


26-043 


AB Aur 


5 


B9.5-A0 


128309 


Ml 


1404 


2 


0.98 


1.00 


954 


2 


0.98 


1.00 


451 1 




0.99 


26-050 


2MJ04554757/801 


2 


M4.75/5.6 


126544 


Ml 


209 


1 




0.74 


118 


1 




0.67 


96 1 




0.90 


26-067 


SU Aur 


2 


G2 


128035 


Ml 


26452 


12 


0.98 


1.00 


6481 


5 


0.99 


1.00 


19971 11 


0.98 


1.00 


26-072 


HBC427 


3 


K7 


128090 


Ml 


13791 


8 


0.98 


1.00 


6147 


6 


0.98 


1.00 


7643 8 


0.97 


1.00 
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Fig. 2. EPIC/pn time-series for three TMC members, demon- 
strating the various shapes of the lightcurves: a 'canonical' flare 
(XEST- 14-006), a smooth decrease of the count rate (XEST-15- 
040), and constant signal (XEST-22-021). The blocks resulting 
from the MLB analysis are overlaid on the binned lightcurve 
(1000 s bins). Data is background subtracted. Time-intervals re- 
moved by the background filter are marked with grey-shades. 
For clarity the background lightcurve is not shown. 



We recall that the parameter describing the variability ac- 
cording to the MLB method (CL) is a threshold, indicating a 
lower limit for the significance (or confidence) with which the 
detected variability is physical, i.e. not spurious. The param- 
eter describing the variability according to the KS test (P) is 
a probability for source variability. Therefore, variability de- 
tected with a low value for CL is compatible with a high value 
of P. 

In general, we find good agreement between the detection 
of variability with the KS test and with the MLB method. In this 
article we distinguish variable from non-variable sources based 
on the results from the MLB analysis for the broad band. There 
are 67 sources with a number of blocks A^b > 1 • For most of 
these sources the KS test probability for variability is very high. 
Exceptions are XEST-19-049 and XEST-08-058, with Pks < 
0.95. After visual inspection of their lightcurves we assigned 
these two sources to the non-variable sample. Then, the total 
number of variable TMC sources is 65, representing 65/126 = 
52 % of all TMC members detected in the XEST. 

In Fig.|2we display the lightcurves of three TMC members, 
representing some of the typical observed shapes: fast 'impul- 
sive' variations on time-scale of hours carrying the signature of 
flares (e.g. XEST- 14-006), gradual variations possibly repre- 



senting a fraction of a long-duration event (e.g. XEST-15-040), 
and constant emission (e.g. XEST-22-021). 

3. Quantifying variability: Definition of flares 

To classify the varia bility we examined the amplitudes and the 
timescales involved. IWolk et aTl (l2005t) have divided the seg- 
ments resulting from a similar MLB analysis in three differ- 
ent types, representing 'characteristic', 'elevated', and 'very 
elevated' intensity. The relatively short exposure time of the 
XEST observations together with the typical strong variability 
of young stars makes it difficult to tell what is the characteristic 
emission level of a given source. Therefore, we follow a some- 
what simplified scheme, in which we consider the segment with 
the lowest count rate in a given time series as the 'characteris- 
tic' emission. The characteristic count rate is henceforth de- 
noted Rch- 

Similar to lWoIk et al. 1 ll2005l) . for the definition of flares we 
make use of the amplitude and the derivative of the segmented 
lightcurve. The amplitude is defined as 
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and the derivative as 

Ai+, = rrrrrr r [ct S ^] 



MIN[ti+i,ti] 



(2) 



In Eq.[2and Eq.|2^i is the count rate of segment /, cr; its un- 
certainty, and ti its duration. 

We define as a flare those variabiHty detections for which 
both of the following two criteria are fullfilled: (i) the ampli- 
tude Ai > 1.5 in one or more consecutive segments, and (i) 
the maximum of the derivatives in these segments exceeds the 
flireshold of A; > 5 10 ^ 

3.1. Results from flare detection 

With the procedure described above flares are detected on 30 
of the 65 variable sources. Two sources have shown more than 
one event (XEST-23-047 wifli two flares, and XEST-26-067 
with three flares). The results of the flare detection process are 
summarized in Table |2 For all flaring sources the XEST ID 
(col. 1), the amplitude of the flare Ap (col. 2) and the maximum 
of the derivative [MAX(A)]p during the flare (col. 3) are given. 
Furthermore, we list the duration rp (col. 4), obtained by sum- 
ming the length of all segments that define the flare. Owing to 
the short exposure times, only 18 events are observed in their 
entirety. For the remaining flares the duration given in Table |2 
is a lower limit. In col. 5 the quiescent luminosity Lch is given, 
and cols. 6 and 7 represent the flare luminosity Lp and flare 
energy Ep. 

Before the conversion to luminosities, the count rates were 
multiplied with a correction factor that accounts for the source 
photons outside the extraction area and for the vignetting. Then 
PIMMS' has been used to obtain the unabsorbed flux for a 
given co unt rate under th e assumption of a 1-T Raymond-Smith 
model fcRavmond & Smiths. 1977 ) subject to photo-absorption. 
For the column density and temperature of the model we 
ad opted the re sults from the spectral fitting presented in Table 6 
of iGiidel et al.. (.2006.) for each individual X-ray source. In the 
case of 2-T spectral fits we computed the emission measure 
weighted mean temperature for the input to PIMMS. The flux 
was converted to luminosity assuming a distance of 140 pc for 
all sources. 

This procedure applied to the characteristic rate R^h yields 
L^-h, while from the average count of all segments that define 
the flare we obtain the average source luminosity during the 
flare. The flare luminosity Lp was obtained by subtracting Lch 
from the average source luminosity during the flare. Finally, the 
flare energy was computed by multiplying Lp with rp. In cases 
where only the rise or only the decay were observed. Zip is a 
lower limit to the energy emitted during the event. 

The last column of Table [^provides a flag that character- 
izes the shape of the lightcurve for the brightest sources with 
flares, adopted from iFranciosini et al. I (I2OO6,) . They present a 
detailed study of the time-evolution of spectral parameters for 
the brightest variable XEST sources, and classify them into 
four groups: The first group is given by typical or 'impulsive' 



' The Portable Interactive Multi-Mission Simulator (PIMMS) is ac- 
cessible at_http://asc.harvard.edu/toolkit/pimms.jsp_ 



flares with short rise and longer decay. Some flares are termed 
'atypical' because they show a gradual rise. The group labeled 
'smooth' in Table |2 is defined by lightcurves that are either 
slowly decaying or slowly rising throughout the whole observa- 
tion. For these sources the characteristic level determined with 
the MLB procedure is most likely too high. We did not attempt 
to classify the flares on the fainter stars. In many cases they are 
described by a single elevated block in the MLB analysis, such 
that their shape is difficult to determine. 



4. Variability statistics 

Table |5] summarizes the variability statistics for different 
classes of young stars: protostars, cTTS, wTTS, brown dwarfs 
(BDs), Herbig Ae/Be (HAeBe) stars, and objects with uncer- 
tain classification. The integer numbers for the object typ e 
identifiers (col. 2) have been introduced bv'Oude l et al.lll2006l) . 
We examined the variability in the soft and hard band, analo- 
gously to the broad band. Note that in some of the subgroups 
the sample size (col. labeled 'Number studied') is slightly 
smaller in the restricted energy bands due to insufficient statis- 
tics. 

In the broad band, variability is detected in about half of 
the time series. This result seems not to depend on the ob- 
ject type. In particular, roughly the same fraction of cTTS and 
wTTS are variable. For the HAeBe stars the number of objects 
is too small to draw statistically valid conclusions. The absence 
of variability on 7 of the 8 detected BDs may be explained by 
their faintness that makes it difficult to discern small-amplitude 
variations (see Sect. 14. 1> . 

The number and the number fraction of variable sources in 
the hard band are similar to those in the broad band, but signif- 
icantly less variations are detected in the soft band (< 30 % of 
the total sample are variable at 0.3 - 1.0 keV). This suggests a 
tight relation between variability and heating processes. On the 
other hand, in some cases the detection of variability in the soft 
band may also be impeded by extinction. The visual absorption 
Ay is on average larger for the sources that show variability in 
the hard but not in the soft band ('group H'), with respect to the 
sources that are variable at both soft and hard energies ('group 
H-hS'). In particular, ~ 40% of the sources from group 'H' 
have Ay > 4 mag, compared to only ~ 10% of the sources 
from group 'Hh-S'. 

Table |4] summarizes the flare statistics for the different 
types of young stars. The flare detection process was run only 
for the broad band. A fraction of 20 - 30 % of each YSO 
class with sufficient statistics has shown detectable flares. The 
flare frequency is marginally larger for cTTS with respect to 
wTTS, consistent w ith earlier results based on ROSAT data 
dStelzer et al. 1120001) . 

Analogous to the diff'erent YSO types, we compare the vari- 
ability and flare statistics of stars with diff'erent spectral types 
in Table |5] Only the results for the broad band are shown. 
The fraction of variable stars drops along the spectral type se- 
quence. However, below we show that this is due to an obser- 
vational bias rather than a physical effect. 
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Table 2. Parameters of X-ray flares detected on TMC members during the XEST: amplitude Ap and derivative [MAX(A)]p are 
used to identify flares; the duration Tp, characteristic luminosity L^h and flare luminosity Lp, and the fl are energy Ep are derive d 
as described in Sect. 13. II the last column provides a flag that classifies the shape of the lightcurve (see iFranciosini et al. LI2006I) . 



XEST 
No. 


Af [MAX(A)]p 
[IQ-^ cts 


Tf 
[ks] 


log Lch 
[erg/s] 


logLF 

[erg/s] 


log^F 

[erg] 


Class 


02-022 


3.0 


21.5 


9.3 


29.3 


29.5 


33.5 




03-022 


2.5 


9.1 


13.5 


29.8 


30.0 


34.1 




04-016 


2.8 


27.0 


> 11.5 


30.7 


30.8 


> 34.8 


smooth 


05-024 


4.8 


90.5 


> 1.5 


30.3 


30.9 


> 34.1 




07-041 


3.7 


98.8 


2.7 


29.7 


30.1 


33.5 




09-026 


1.9 


6.0 


> 4.4 


31.0 


31.0 


> 34.6 


smooth 


11-057 


6.3 


48.2 


> 16.6 


30.4 


30.9 


> 35.1 


smooth 


13-004 


3.4 


45.0 


> 7.2 


29.3 


29.6 


> 33.5 




14-006 


5.7 


123.0 


15.0 


29.8 


30.2 


34.4 




15-040 


2.6 


5.8 


> 20.7 


30.8 


30.8 


> 35.2 


smooth 


17-058 


2.6 


7.2 


> 5.5 


29.6 


29.7 


> 33.5 




17-066 


3.0 


95.5 


9.5 


30.3 


30.6 


34.6 


atypical 


18-004 


4.0 


31.1 


> 6.0 


29.7 


30.2 


> 34.0 




20-005 


1.8 


7.1 


> 1.2 


30.7 


30.6 


> 33.7 




20-022 


80.3 


4121.5 


3.1 


30.4 


32.1 


35.6 




20-056 


3.6 


58.4 


> 7.2 


30.1 


30.5 


> 34.3 




22-047 


4.0 


12.7 


> 41.6 


30.1 


30.4 


> 35.1 


smooth 


22-089 


2.1 


9.0 


> 4.2 


30.3 


30.3 


> 34.0 


smooth 


23-002/24-002 


5.6 


69.3 


15.2 


29.4 


30.0 


34.2 




23-032/24-028 


2.9 


39.5 


46.6 


30.7 


30.7 


35.4 


atypical 


23-033/24-029 


8.7 


64.9 


8.6 


29.2 


30.0 


33.9 




23-045/24-038 


2.1 


150.5 


> 26.9 


29.6 


29.6 


> 34.0 




23-047/24-040 


1.6 


40.3 


40.9 


30.9 


30.7 


35.3 


atypical 


23-047/24-040 


8.2 


174.6 


> 17.5 


30.9 


31.6 


> 35.9 


atypical 


23-048 


7.1 


31.3 


> 2.7 


28.7 


29.6 


> 33.0 




23-050/24-042 


13.0 


637.8 


47.4 


30.0 


30.6 


35.3 


impulsive 


23-063/24-055 


3.1 


11.1 


36.0 


29.7 


29.9 


34.5 




23-074/24-061 


1.8 


8.5 


6.7 


30.4 


30.3 


34.2 


impulsive 


26-067 


3.0 


35.8 


19.9 


31.1 


31.3 


35.6 




26-067 


2.0 


23.9 


9.0 


31.1 


31.0 


35.0 




26-067 


1.7 


18.3 


13.0 


31.1 


30.9 


35.0 




26-072 


5.2 


87.6 


56.0 


30.6 


30.9 


35.6 


impulsive 


28-100 


3.2 


62.2 


23.8 


30.2 


30.4 


34.8 


impulsive 



Table 4. XEST X-ray flare statistics for different types of 
young stars in the 0.3 - 7.8 keV broad band derived with the 
MLB technique. For each YSO class the number of detected 
sources and the number of flares are given, and the fraction of 
flares with respect to the total sample. 



Object 


Type 


— Members — 


Fraction 


type 


identifier 


detected 


flare 


flare 


Protostar 


1 


9 


2 


22% 


cTTS 


2 


52 


16 


31% 


wTTS 


3 


51 


11 


22% 


Brown dwarf 


4 


8 





0% 


HAeBe 


5 


2 


2 


100% 


Unknown 


9 


4 


2 


50% 


Total 




126 


33 


26% 



Table 5. XEST X-ray variability statistics for diff'erent spectral 
types in the 0.3 - 7.8 keV broad band; see Tables|3]and|4]and 
text in Sect.|4]for details. 

Spectral Members Fraction 



type 


detected 


variable 


flare 


variable 


flare 


BAF 


2 


2 


2 


100% 


100% 


G 


4 


4 


3 


100% 


75% 


K 


35 


23 


12 


66% 


34% 


M 


78 


32 


14 


41 % 


18% 


Unknown 


7 


4 


2 


57% 


29% 


Total 


126 


65 


33 


52% 


26% 



4.1. Observational biases 

Next to residual uncertainties in the background subtraction, 
the capability to detect variations is affected by Poisson statis- 
tics and the inhomogeneous length of the XEST observations. 
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Table 3. XEST X-ray variability statistics for different types of young stars in the broad, soft and hard band derived with the 
MLB technique. For each energy band the number of detected sources, the number of variable sources, and the number fraction 
of variables are given; see text for details. 







Broad band: 0.3 - 


7.8 keV 


Soft band: 0.3 - 


1.0 keV 


Hard band: 1.0- 


7.8 keV 


Object 


Type 


Number 


Fraction 


Number 


Fraction 


Number 


Fraction 


type 


identifier 


detected 


variable 


variable 


studied 


variable 


variable 


studied 


variable 


variable 


Protostar 


1 


9 


4 


44% 


7 





0% 


9 


4 


44% 


cTTS 


2 


52 


28 


54% 


50 


9 


18% 


52 


27 


52% 


wTTS 


3 


51 


26 


51% 


51 


19 


37% 


51 


29 


57% 


Brown dwarf 


4 


8 


1 


13% 


7 





0% 


7 





0% 


HAeBe 


5 


2 


2 


100% 


2 


2 


100% 


2 


1 


50% 


Unknown 


9 


4 


4 


100% 


















Total 




126 


65 


52% 


117 


30 


26 % 


121 


61 


50% 



The resulting biases are discussed by investigating the ampli- 
tudes of the blocked lightcurves. 

We define the amplitude of the lightcurves in absolute terms 
as Aabs - RmsLx - Rch, whcrc 7?,nax IS the highest block count 
rate. The left diagram in Fig.|3shows Aabs of all variable stars 
with known spectral type as a function of the characteristic rate. 
Clearly, the absolute amplitude is correlated with Rch- The ex- 
istence of a lower envelope is due to the statistics that impose 
a sensitivity limit on the detection of variability. In order to be 
recognized as variable, the signal in a given block must ex- 
ceed a certain signal-to-noise level, k: 



XEST-20-022 



ch.i 



(3) 



where A^i is the observed number of counts in block /, A^^h.i is 
the number of counts in block / that represent the characteris- 
tic level, and cr is the standard deviation as a measure for the 
noise. Transformed to count rates, this yields for the absolute 
amplitude 



1 



ch,max) ^ ^ \l 

^max 



(4) 



0.14 
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0.10 
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0.04 
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This sensitivity limit is approximated by the dashed line in 
Fig. Its location is not uniquely determined, because it de- 
pends on the duration of the block with the maximum count 
rate. For reasonable values of fmax ~ 5...20ksec the plotted line 
corresponds to k ~ 3.5. ...7. 

On the right hand side of Fig. |3] the relative amplitude, 
Arei = (^max - Rch) / Rch IS shown. In this representation the 
sensitivity threshold has a functional dependence of 1 / V^ch- 
Therefore, events that are small in relative terms are more eas- 
ily detected on bright stars, while events that are small in abso- 
lute terms are more easily detected on faint stars. 

An approximate empirical upper envelope to the observed 
variations is marked in Fig. |3] with a dotted line. In contrast to 
the lower threshold that is - as explained above - induced by 
our sensitivity limit, the absence of large amplitudes is prob- 
ably the result of the intrinsically rare occurrence of strong 
intensity changes combined with our limited observing time. 
During the XEST the majority of TMC members did not un- 
dergo variations by more than a factor 8 in count rate. 

The different plotting symbols in the right diagram of 
Fig. 13 distinguish not spectral types but the YSO classes 
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Fig. 4. Background subtracted broad band lightcurve of the 
protostar IRAS-04108-H2803B (XEST-20-022) showing a 
large flare. The data is dominated by the hard band. The soft 
emission is absorbed due to high column density. 



from Table |3 An exceptionally large amplitude is shown by 
IRAS 04 108 +2803 B (XEST-20-022), a Class I protostar This 
object has shown a spectacular flare near the end of the observa- 
tion, with a rise in count rate by a factor of ~ 80. Its broad band 
lightcurve is displayed in Fig.0] The pre-flare quiescent count 
rateofIRAS04108-H2803B is - 1.9 lO^^ cps in EPIC/pn, near 
the detection limit of th e observation . Accor ding to the spectral 
analysis carried out by Giide l et'aP ( l2006l) the source suffers 
from strong extinction (log A^H[cm~^] ~ 22.9). Indeed, the soft 
emission is completely absorbed, explaining the sharp contrast. 

Combining the statistical sensitivity threshold with the em- 
pirical maximum of the observed amplitudes, the dynamical 
range for detectable variations is much larger for bright stars 
than for faint stars. Therefore, one expects to find more vari- 
ables among stars with high characteristic emission. To demon- 
strate the distribution of R^h in the total sample, at the bottom 
of Fig.|3the positions of non-variable stars of the different sub- 
groups are indicated. For Ret > 0.3 cts ' about 50% of the 
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Fig. 3. Absolute amplitude (left) and relative amplitude (right) of the segmented lightcurve vs. characteristic count rate for 
variable stars. On the left the different plotting symbols distinguish stars with different spectral type, on the right stars of the 
different YSO classes. The dashed line denotes the sensitivity limit according to Poisson statistics for k / VWx - 0.05. The 
dotted line marks the approximate location of an empirical upper envelope for the amplitudes at /?n,ax = 8 /?ch- Upward pointing 
arrows are positioned at the average count rate of non-variable stars of the different spectral type and YSO groups, respectively. 



sample is variable. Similar observations have been discussed 
in the literature, e.g. Stelzer etal. (2005), JFlaccomio et al. 

have shown that for a sample in the field of a given X- 
ray observation all stars above a certain number of counts are 
variable. Probably longer observations would yield the same 
result for the TMC sample. 

4.2. Variability on different groups of pre-MS stars 

Taking account of the biases discussed above, we can proceed 
to a further examination of the variability on different subsam- 
ples of TMC members. A comparison of Table |5] with Fig. |3 
(left) suggests that the smaller fraction of variables among M 
stars when compared to K stars is a result of their fainter char- 
acteristic emission. Indeed, the vast majority of M stars gather 
at /?ch < 0.01 cts ' , where the detection of small variations is 
impeded by Poisson statistics. The K stars follow a wider distri- 
bution including higher values of Rch, such that the detection of 
variability is favored. When only sources with R^h < 0.02 ct s ' 
are considered, the fraction of variable sources is 46 + 19 % for 
the K stars, and the fraction of variables is 50 + 12 % for the 
M stars, i.e. the variability statistics for the two spectral type 
classes are indistinguishable. 

In an analogous way, the variability statistics of cTTS and 
wTTS c an be compared. It was repeatedly shown in t he lit- 
erature ( iNeuhauser et al. L Il995[ IStelzer & Neuhauser L l200ll) 
that wTTS are on average X-ray brighte r than CTTS. Th is has 
also been confirmed for the XEST (.Giidel et al.L l2006l) . and 



therefore one may expect to find a reduced fraction of variable 
sources among cTTS with respect to wTTS. However, from 
Table |3 it results that a similar fraction of cTTS and wTTS 
(~ 50 %) is variable. As seen from Fig. (right) both sub- 
groups have a large fraction of faint objects on which vari- 
ability is diflicult to detect (~ 65 % of all cTTS and 45 % of 
all wTTS have R^h < 0.02 cts '). In the faint subsample with 
Rch < 0.02 cts"' the fraction of variable sources among the 
cTTS is 41 ± 11 %, and the fraction of variables among the 
wTTS is 26 + 1 1 %. This is a significant reduction with re- 
spect to the sample without restriction in Rch, but the results 
for cTTS and wTTS are undistinguishable within the statistical 
uncertainties. 



From Fig. |3it can also be seen that the non-detection of 
variability on most of the BDs is likely to be related to their low 
characteristic count rates: The XEST is s ensitive only to large , 
and presumably rare, flares on the BDs. lOrosso et al. I (l2006l) 
have identified a flare in the binned lightcurve of CFHT-BD- 
Taul (XEST-17-068). Our MLB algorithm does not recover 
this event, although there is some evidence for variability in 
the binned lightcurve. According to the KS test the broad band 
photon time series of XEST-17-068 is variable at the 99 % con- 
fidence level. Possibly this flare was missed in the MLB anal- 
ysis due to high background. Its broad band lightcurve, dis- 
played in Fig.|5l shows that part of the flare was removed from 
the data with our high-background filter. 
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Fig. 5. Background subtracted broad band EPlC/pn lightcurve 
of the BD CFHT-BD-Tau 1 ( = XEST-17-068) . There is some 
evidence for the possible flare reported bv .Grosso et al. 
at f~ 7000- 10000 s. 



5. Flare duration and frequency 

The duration of the observed flares ranges from a few ksec 
to ~ 60 ksec, with a clustering of events near ~ 10 ksec. 
Obviously, the XEST is biased against the detection of long- 
duration events, due to the short exposure times (typically 
~ 30 ksec). Recall that the distribution of flare durations de- 
rived from the 13d-long COUP exposure JWolk et al. L l2005l) 
ranges from 1 hr to 3 d with a peak near 18 h. Consequently - 
leaving a part subtle d i fferenc es in the definition of a flare in the 
work bv IWolk et al. I (|2005) and the one presented here, and 
ignoring that the two studies concern stars of different mass 
range that may influence the flare characteristics - it is no sur- 
prise that nearly half of the flares we observed in the XEST are 
only partly within the observing time. As mentioned above we 
can give only lower limits on the duration of these events. 

The frequency of flares for the total XEST sample of TMC 
members computed from the total observing time T^bs and the 
number of observed flares is 1 flare per star in 200 ksec. In com- 
puting this number we have assumed that we are sensitive to the 
detection of flares throughout the whole observing time. This 
should be a reasonable approximation, because observing gaps 
due to high background are mostly shorter than the typical du- 
ration of the flares. A serious limitation, however, comes from 
the fact that we are comparing stars whose characteristic count 
rates differ by as much as a factor of 1000. As a consequence 
the events detected in our observations form an incomplete sub- 
sample of the total flare population. In Sect. |6l we argue that 
the flare detection is complete for events with Ep > 10^^ erg/s. 
If only those events are considered the flare frequency of the 
TMC members reduces to 1 flare in 770 ksec. Our flare detec- 
tion is also biased by the limited length of the XEST obser- 
vations, typically ~ 30 ksec. Indeed, during the long (115 ksec) 
merged observation XEST-23/XEST-24 8 flares were identified 
on the 18 stars detected in both exposures XEST-23 and XEST- 
24. Counting only the events with energy above lO-'^erg the 
XEST-23/XEST-24 field has a flare frequency of 1 in 520 ksec. 



These numbers can be compared to the flare frequency 
observed during the COUP for the young solar-analogs in 
the ONC. There are 30 flares that emit more than 10^^ erg/s 
(Table 6 of Wolk etaLl l2005h . and - following the discussion 
in the COUP article - we estimate that their occurrence rate is 
~ 1 /900ksec, similar to our result for a sample of wider mass 
range in the TMC. 

The fraction of the total observing time during which a 
given group of stars was found in the flare state is defined by 



(5) 



Tobs 



where the sum goes over the durations of all flares. For the to- 
tal TMC sample we find Fp ~ For the subgroups of cTTS 
and wTTS the fractional time in the flare state is found indis- 
tinguishable in both cases at ~ 8-9%. IStelzer et al. I feOOOh 
have found ~ 1 % for the flare rate of TMC stars, evaluating 
pointed ROSAT observations. This discrepancy may indicate 
that ROSAT missed a substantial number of flares due to gaps 
in the observing sequence resulting from Earth occultations, or 
due to uncertainties in the flare durations related to these gaps, 
or due to the smaller effective area or missing sensitivity for 
more energetic photons. 

6. Flares and the coronal heating process 

To study the effect of the flare population on coronal heating, 
we now construct the distribution of flare occurrence rate in 
total released energy. For the Su n, it was shown that this distri- 
bution foUows a power law (e . g. iHudson L 1 1 9911) 



dN 

Ie 



(6) 



The radiated energy released in flares in the observed energy 
band is obtained by integrating this differential distribution. 
Generally, observational studies are affected by a sensitivity 
limit that impedes the detection of very small flares. It is com- 
mon practice to extrapolate the power-law to energies below the 
detection threshold in order to get hold of the energy in small 
unresolved events. Such 'nano-flares' are assumed to exist be- 
cause they can be considered as the heating agent that gives 
rise to the quiescent corona of the Sun and other magnetically 
active stars. If the power-law index a > 2, very small flares can 
in principle contribute an unlimited amount of energy, because 
the integral of Eq. |6l diverges for /imin — > 0. Clearly, in such 
cases there must either be a lower cut-off to the flare energies, 
or the power law must turn over to a < 2 at low energies. 

Although we are limited to measurements exceeding 
log /ip [erg] « 33.5 in our survey, we note that a power- 
law distribution has been followed down to luminosities of 
log Ep [erg] = 30.5 for a few examples of somewhat more 
evolved, magnetically active low-mass stars and to levels 
as low as log Ep [erg] - 25 in the solar corona. For the 
Sun the situation has remained inconclu sive, with a between 
1.6 — 2.6 over a vvide range of energie s (Crosbv etal. I ll 9931: 
iKrucker & Benz L Il998t lAschwandenet al. . .2000 ). For stars 
the hard X-ray range, mostly used in solar flare observations. 
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is inaccessible. In recent statistical studies of EUV flare en- 
ergy d istri butions on magnetica lly active stars lAudard et al. I 
(1 1999*) and' Audard et al. I (I2OOO') find power-law indices close 
to the critical value (a ~ 2). Further observ ations of ste llar flare 
ene rgy distributions are summarized bv lKashvap et alTI (l2002l) 
and lOiidel et al. I ll2003l) . who point out that the slope may de- 
pend on the spectral range of the instrument. A comparison 
of the distributions observed in the EUV, in soft, and in hard 
X-rays showed that these diff'erent temperature regimes corre- 
spond to different energy ranges from 'nano'- to 'milli-flares', 
and that they represe nt different density, press ure, and emis- 
sion measure regimes JAschwanden etal.Ll2000 ). Therefore, it 
is unclear if the extrapolation from the large observable flares 
towards the lowest energies is valid. 

The analysis of such distributions is subject to several 
caveats. One of them is that small flares overlapping larger 
ones are masked a nd missed in the flare detection process. 
lAudard et al. llEoOO) have mediated this problem by applying a 
correction factor to the flare rate at each energy. This correction 
takes account of the fact that the 'effective' time available for 
the identification of small flares is given by the total observing 
time (Tobs) reduced by the sum of the durations of larger flares 
(t>e), such that the flare rate at energy E must be multiplied by 
a factor 
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Fig. 6. Cumulative distribution of flare energies for the pre-MS 
sample in the TMC observed during the XEST (Table |2j and 
for the sola r-analogs of the O NC observed during the COUP 
(Table 6 of lWolket ail 12005*). The dotted lines represent the 
ML estimates for the slope of the distribution log = j6-log £-1- 
C in the range of energies above the cutoff value E^m indicated 
by the open circle. 



/e = Tohs/iTobs, - t">e). (7) 

6.1. Flare energy distributions: TMC and ONC 

Systematic studies of the energy release during flares on pre- 
MS stars have been carried out for only two star forming re- 
gions so far: the TMC (this study) and the solar-analogs of 
the ONC dWolk et al. L l2005h . In Fig.|6lthe cumulative num- 
ber distribution [A'(> Ep)] of flares observed in both samples 
are shown in double logarithmic form. 

We derived the distributi on of the ON C from the mean flare 
energies given in Table 6 of lWolk et al.1 (12005). In contrast to 
our XEST study of the TMC, these energies were computed 
from the spectrum during the flare state. In some cases the val- 
ues for Ep are dubious as a result of poor spectral fits, and 
we exclu ded events flagged with 'c', 'e', or 'f in Table 6 of 
IWolk et al. . (,2005 ) from the analysis of the energy distribution. 

We performed for both distributions the correction for over- 
lapping flares described above by multiplying each value of 
A^(> E) with its appropriate factor /e. For the ONC data we 
assumed a total observing time of 850 ksec ■ 22, where 22 is the 
number of stars in the sample after removing the 5 stars with 
poorly defined flares mentioned above. The duration of each 
flare is taken from Table 6 of Wolk et al. ( 2005). The effect of 
the correction factor /e is to steepen the slope, but in practice 
its influence is negligible because the flare frequency is low, 
such that T>E is only a minor fraction of Tabs- 

The faintest flare detected in the ONC sample has an en- 
ergy of 10^"^^ erg/s, while in the TMC we have access to events 
that are one order of magnitude lower. This difference in sen- 
sitivity corresponds roughly to the difference in X-ray flux for 
a source of given luminosity that arises from the distance ratio 
(</oNc/i/TMc)^ (140 pc for TMC vs. 450pc for the ONC). 



6.2. Determination of the power-law slope 

As discussed above, the distribution of flare energies is ex- 
pected to follow a power-law. However, the observed distribu- 
tions flatten towards lower energies. This saturation is prob- 
ably an observational bias, indicating that smaller events are 
missed by the flare detection process. On the other hand, the 
steep high-energy portion represents large flares that are read- 
ily detectable, such that above a certain threshold energy E^ut 
the flare identification is likely complete. 

We have determined both the slope /3 of the function 
log A^(> Ep) - -/3 ■ log Ep + C and the cutoff energy Ect with 
a maximum likelihood method. According to ICrawford et al. I 
( 1970) the ML estimate for the power law slope is 

l = iyin(A) (8) 

where is the total number of flares above E^m, and are 
their energies. To test if Eq. |S]is a good representation of the 
data we have performed the transformation 

y, = (l-(A)-^)/(l-(f^)-'') (9) 

^cut ^cut 

and tested this distribution a gainst departure from a uniform 
distribution using the KS test dCrawford et al. [[19701) . 

We have evaluated Eq. |8] and |9l for a range of cutoff en- 
ergies. The result is shown in Fig. It presents in the same 
graph, but with different scales, the run of /3 (filled circles) 
and the probability Pks from the KS test (open squares). High 
values of f ks indicate that the distribution of y-, is compatible 
with the null hypothesis of being uniform in the interval [0,1]. 
Therefore, the power law approximation of the observed flare 
energy distribution is justified for cutoff energies in the peak 
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Fig. 7. Result from the ML estimate of the power-law index for 
the flare energy distribution of the XEST sample in the TMC 
and the COUP sample in the ONC as a function of minimum 
flare energy included {Ecm)- filled circles - slope /3 of the double 
logarithmic cumulative distribution of flare energies (log = 
-/3 ■ \ogE + C); open squares - probability of KS test that the 
above power law is a good approximation of the observed flare 
energy distribution. 



TMC sample, e.g. the definition of the flares and the evaluation 
of the flare energies. In our analysis we have ignored that the 
observed energies for some flares in the XEST are lower limits. 
No information on a possible incomplete coverage of the ONC 
flares is given bv lWolket al.lhOQSl) . However, as judging from 
their Fig. 3, some of the events should be flagged lower limits 
because data gaps during the flare have caused a loss of pho- 
tons for the spectral analysis. We have verified in simulations 
that the power-law index /? is not affected by ignoring the cen- 
soring of some energies. If a certain number of values at ran- 
dom drawn from a power-law distribution are decreased by a 
random fraction of their original value, this results in a left- 
ward shift of the cumulative number distribution, but the slope 
remains unmodified. This means that we implicitly assume a 
random distribution of the lower limits across the range of ob- 
served flare energies. While in practice this assumption remains 
unproven, this is unlikely to be of major importance to our re- 
sults. 

We also note that we did not perform a thermal modeling 
of individual flares to estimate their radiated energies in detail, 
given that the spectra were too faint in many cases to obtain 
reliable temperature, emission measure, and absorption param- 
eters. It was implicitly assumed that the flare parameters are 
similar to the coronal parameters determined for the overall 
emission. Given the relatively moderate dynamic range of flare 
amplitudes and also their moderate amplitudes compared with 
the overall emission level in the respective light curves, we be- 
lieve that this approximation is appropriate. 



6.4. Completeness limit for flare detection 



of Pks (Ecut)- From Fig. |51the power law index /5 is expected 
to increase with increasing Ecui- For cutoff energies where the 
flare energy distribution is well represented by a power law (i.e. 
for high values of Pks), the flare detection process is prob- 
ably complete. Therefore, further increasing Ecut should not 
change the slope, and /3 (Ecut) reaches a plateau. Obviously, at 
the largest cutoff energies the result is compromised by low- 
number statistics. The number of flares with Ep > E^m are 
given on the top of both panels in Fig.0 

6.3. Results and caveats 

The best match for the power-law index determined from Fig.0 
is jStmc = 1.4 + 0.5 for the TMC and /3onc = 0.9 ± 0.2 for the 
ONC. These slopes are plotted in Fig.|6l and the corresponding 
EcM is marked with an open circle. It suggests that the flare 
detection is complete for events of log Ep [erg] > 34.9 in the 
case of the TMC, and logEp [erg] > 35.3 for the ONC. The 
cumulative distribution of Fig. |S] has the derivative dN/dEp ~ 
E^", where a - \ + [5. Therefore, from our observations we 
find QfTMC ~ 2.4 and aoNC ~ 1-9. Both values are consistent 
with each other and with a > 2, considering the substantial 
uncertainties. 

Besides the statistical errors, there are various uncertain- 
ties and differences in the approach between the ONC and the 



We estimate the minimum signal for the detection of a 
flare with edge energy fip.edge = 10^^ erg and a typical du- 
ration of lOksec. For such events the flare luminosity is 
logLF,edge[erg/s] - 31. We compute the value of Lp.edge/^'ch 
that corresponds to this threshold from the characteristic count 
rate of each star. Agai n, we use PI MMS with spectral param- 
eters from Table 6 of iGiidel et alJ ( 1200 6) to determine an in- 
dividual count-to-flux conversion factor for each source, that 
allows to determine L(,h- 

The resulting amplitude Lp.edge/^'ch is visualized in Fig.|8l 
The choice of upward pointing arrows as plotting symbol is 
motivated by the indications discussed in Sect. 16.21 that all 
flares with energy above the cutoff fip.edge are detected, such 
that Lp.edge/i'ch represents a completeness threshold. As can be 
seen, for given R^h the differences in the spectral shape result in 
a spread of more than one order of magnitude in the luminos- 
ity amplitude. This is illustrated by the dotted lines that repre- 
sent two different conversion factors. For comparison in Fig.|8] 
the amplitudes of all observed flares are also shown (asterisks). 
Generally, the 'completeness limit' corresponds to much higher 
amplitudes than observed, leaving space for a large popula- 
tion of undetected events. Filled asterisks denote those flares 
with Ep > 10^^ erg. As expected most of these are well above 
the derived completeness limit for their host star (marked with 
thicker arrows). In some cases they are below, and this can 
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Fig. 8. Flare sensitivity limit. Shown is for each star the am- 
plitude above which all flares with an assumed emitted energy 
of 10^''' erg/s and a duration of lOksec are detectable (upward 
pointing arrows). The conversion from count rate to luminosity 
depends on the spectral shape, and is expressed by the dotted 
lines that indicate two different count-to-energy conversion fac- 
tors. Asterisks denote the amplitudes of observed flares; filled 
asterisks are those flares with E > 10^^ erg/s, and the flare sen- 
sitivity limit of the corresponding star is shown as thicker ar- 
row. 



be explained by their duration being much longer than our as- 
sumption of lOksec. 

7. Variability of cTTS and wTTS 

A comparison of the variability characteristics of cTTS and 
wTTS is of interest because of the possibly different emission 
mechanisms: Next to the canonical interpretation of magnetic 
structures with footpoints anchored on the star, cTTS may pro- 
duce X-ray emission related to the accretion process and/or in 
magnetic structures connecting the star with the circumstellar 
disk. In wTTS the absence of a disk and accreting material im- 
ply that the most plausible explanation of the short-term X-ray 
variability is magnetic flaring. 

Observationally it is not straightforward to distinguish be- 
tween variability originating from these different causes, unless 
there is sufficient signal for detailed modelling of the spec- 
tral evolution that gives access to the strength of the mag- 
netic field an d the length scale of magnetic structures (see e.g. 
iFavata et al. L 12005 ). A few of the brightest XEST sources are 
examined in this respect by Franciosini et al. ( 2006). Here we 
compare the variability of cTTS and wTTS using a statistical 
approach. 

In Sect. 12 we have established that a similar fraction of 
cTTS and wTTS are variable in the broad band. In the soft 
band, variability is more frequent on wTTS than on cTTS: 
18 + 6% variables among cTTS vs. 37 + 9 % variable wTTS. 
A possible cause for this difference is the on average larger 
extinction of the former ones that removes soft photons from 
the spectrum. About 30 % of all cTTS have Ay < 1 mag, com- 
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Fig. 9. Flare energy and flare duration for cTTS and wTTS. 
Lower limits for flares that are unconstrained because partially 
outside the observing window are indicated with arrows. 



pared to more than two-thirds of the wTTS sample. This trend 
is even more pronounced among the stars that are constant in 
the soft band, and supports the suspicion that soft absorption 
puts stronger limits on the detection of variability in cTTS than 
in wTTS. On the other hand, if a significant amount of X-rays 
from cTTS were produced from accretion, then at least those 
cTTS that are not absorbed should show more - not less - soft 
variability than wTTS, because the temperatures in a T Tauri 
accretion shock (of a few MK at most) can produce only soft 
photons (< 1 keV). This is not observed. An alternative expla- 
nation would be that accretion is steady enough not to produce 
X-ray variability, or that variations in the density of the accret- 
ing material are smoo thed out during the X-ray generation pro- 
cess ( lLamzinL [l 99^ . 

The XEST observations show further that flares (character- 
ized by large amplitude and fast rise) are identified marginally 
more frequently on cTTS than on wTTS (3 1 + 7 % vs 22 ± 7 %). 
A similar conclusion, but likewise not on a strong statistical 
basis, was obtained by Stelzer et al. (2000) based on ROSAT 
observations. In Fig.|9]the flare durations and energies of cTTS 
and wTTS measured during the XEST are shown. There is not 
a clear distinction between the two groups, pointing at a similar 
- most likely coronal - origin for the X-rays from both classes 
of stars. 

8. Conclusions 

We have presented a systematic analysis of X-ray variability 
on pre-MS stars in the TMC based on the XEST. This project 
includes a significant fraction of the known TMC members, 
and presents the deepest X-ray survey so far in the Taurus 
star-forming region. We have used an automated procedure 
for detecting variability and recognizing flares. Flares are de- 
fined with criteria that take into account the amplitude and the 
derivative of the segmented lightcurves. Variability is found in 
roughly half of the sample in the broad band from 0.3-7.8 keV. 
Our analysis in different energy bands has shown that it is im- 
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portant to include hard energies: Variability is detected more 
frequently in the hard band (50 % of sources are variable at 
> 1 keV) than in the soft band (26 % of sources are variable at 
< 1 keV) where, in some sources, strong extinction lowers the 
signal to the point that it becomes difficult to reveal variations. 
Next to this bias there is probably a physical reason behind the 
fact that variability shows up stronger at harder energies, re- 
lated to the notion that variability goes along with additional 
heating. 

The detection of variability is biased because the Poisson 
statistics impose different sensitivity thresholds for sources 
with different count rate, and the sample has a large range of 
count rates (four orders of magnitude). We have estimated the 
influence of this bias on the variability statistics, and evalu- 
ated the variability statistics for different classes of pre-main 
sequence stars (protostars, cTTS, wTTS, brown dwarfs), and 
for different spectral type ranges that roughly correspond to 
mass for cool pre-MS stars. 

Variability is found slightly more frequently on cTTS than 
on wTTS. The difference is only marginal, but when combined 
with similar indications fro m various star-forming re gions dis- 
cussed in the literature (e.g. iFlaccomio et al. LboOfil) . this find- 
ing may deserve further attention. On the other hand, the flare 
energies and durations of cTTS and wTTS are similar, and 
there is no evidence of extra soft variability from accretion 
in cTTS. All in all, the XEST data does not yield support for 
models that predict X-ray emission triggered by reconnection 
relate d to star-d isk magnetospheric interaction and accretion 
(e.g. IShuetal.L ri994: Goodson & Winglee L fl999). although 
such sc enarios may apply in ind ividual cases; see also discus- 
sion bv lFranciosini et alH ( l2006l) . 

There is no dependence of the fraction of var iable sources 
on s pectral type (comparing K and M stars). ICalvet et al. I 
( 1200 4) summarized a correlation between mass accretion rate, 
M, and stellar mass, M, that is found to be M oc M' -^^. If X-ray 
production was related to the accretion activity (e.g., through 
accretion shock heating or through influence on the coronal 
magnetic structure), then stronger accretion would be expected 
to produce stronger (absolute) fluctuations, and therefore X-ray 
variability should be different not only between accreting and 
non-accreting stars, but also between lower- and higher-mass 
cTTS in the sample. Since M stars and K stars show no distinc- 
tion with regard to absolute variability, we find no evidence that 
the accretion pr o cess s ignificantly influences X-ray variability. 
IPreibisch et al. I (l2005l) arrived at similar conclusions for the 
ONC sample based on a comparison of the thermal properties 
ofcTTS and wTTS. 

We have evaluated the flare rates for the events on TMC 
members identified in this work, and compared them to equiv- 
alent data for the ONC sample obtained by Chandra within the 
COUP. In both cases the sample represents pre-MS stars, but 
the ONC study concerns a narrow mass range (0.9 < M < 
1 .2 Mq), while in this work we have included all TMC mem- 
bers detected during the XEST. The frequency of large flares 
(iip > 10^^ erg) in the TMC is 1 event per star in 770ksec, 
roughly comparable to the estimate for the young solar analogs 
in the ONC. 



We have found twofold observational support for coronal 
heating by flares in our stellar sample. First, we have noted that 
the flare amplitude is correlated with the characteristic emis- 
sion level of the flaring star In particular, we found an upper 
envelope to the flare amplitudes as a function of /?ch- While the 
envelope is unlikely to be strict because longer observations 
could detect some of the very rare but very large flares on any 
star, the trend can naturally be explained if flares ar e related 
to the coronal heating process (" Audard et all l2000l) : Either, 
both the flare rate and the overall heating are the result of a 
similar process, probably magnetic reconnection, and therefore 
both the flare rate for a given flare amplitude and the total X- 
ray luminosity are expected to scale with the coronal volume 
available for these processes. Alternatively, the characteristic 
emission could be the direct result of the heating process of a 
large number of unresolved flares occurring continuously on 
the star Both assumptions relate the physical processes rele- 
vant to flares to the overall stellar radiative X-ray losses. In both 
cases, one expects (for a given observing time) that the maxi- 
mum of the absolute amplitude be proportional to R^h, and this 
is indeed what we find (Fig.|3]left). 

Secondly, we have analyzed the cumulative number distri- 
bution of flare energies, making use of our new results for the 
TMC and of the flare energies tabulated bvl Wolket alltosl) 
for the ONC. Our ML estimates show that the distributions of 
both samples can be represented by a power-law, with an index 
ofa- 2.4 ±0.5 for the differential distribution of the TMC, and 
ff = 1.9 ± 0.2 for the ONC. Wolk et al. (2005) had provided a 
value of a ~ 1.7 for the same ONC data, using a linear least- 
squares fit to the double logarithmic representation of the cu- 
mulative number distribution. However, that approach does not 
take into account that the individual bins are not independent. 
Therefore, our method is more rigorous. The power-law slope 
for the number distribution of flare energies in the TMC and in 
the ONC, although not well constrained, is in agreement with 
the range found in previous studies of solar and stellar flares 
(see references in Sect.|6}, indicating that the energy release in 
cool stars is a fairly universal process. The observed values for 
a are near and probably larger than the critical limit (o-crit = 2) . 
Therefore, although the uncertainties preclude definite conclu- 
sions, micro-flare heating may well play an important role in 
heating these pre-MS coronae. 

In summary, the entire sample supports an important role 
of flare-like events in the generation of the X-ray light curves, 
in a similar manner as observed before in non-accreting mag- 
netically active stars. In particular, we conclude that cTTS and 
wTTS reveal similar variability characteristics. These results 
point at magnetic energy release playing the dominant role in 
the generation of X-rays, while there is little evidence for an 
influence of accretion disks (via long-scale magnetic fields) or 
accretion itself (via accretion-shock heating). While neither of 
the latter two can be excluded to contribute to the observed 
X-rays, we claim that neither is responsible for strong X-ray 
losses, but that compact magnetic field annihilation dominates 
the X-ray production. 
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